Full first-principles theory of spin relaxation in group-IV materials 
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We present a generally applicable parameter-free first-principles method to determine electronic 
spin relaxation times and apply it to the technologically important group-IV materials silicon, di- 
amond and graphite. We concentrate on the Elliott- Yafet mechanism, where spin relaxation is 
induced by momentum scattering off phonons and impurities. In silicon, we find a ~ T~ A tem- 
perature dependence of the phonon-limited spin relaxation time Ti and a value of 4.3 ns at room 
temperature, in agreement with experiments. For the phonon-dominated regime in diamond and 
graphite, we predict a stronger ~ T~ 5 and ~ y- 4 - 5 dependence that limits T\ (300 K) to 180 and 
5.8 ns, respectively. A key aspect of this study is that the parameter-free nature of our approach 
provides a method to study the effect of any type of impurity or defect on spin-transport. Further- 
more we find that the spin-mix amplitude in silicon does not follow the Eg 2 band gap dependence 
usually assigned to III-V semiconductors but follows a much weaker and opposite Eg 67 dependence. 
This dependence should be taken into account when constructing silicon spin transport models. 



The physical roadblocks looming in the charge-based 
semiconductor device technology require paradigm- 
shifting approaches to create new logic devices capable of 
lower power consumption and higher performance. This 
has motivated a search for new alternative logic variables, 
among which the spin of electrons is a natural candi- 
date, which needs to be efficiently and reliably injected, 
transported, and detected. Although extensive studies 
have been done in direct-gap materials, the understand- 
ing of spin-transport and spin-life time dependence in 
the technologically relevant group-IV materials is sur- 
prisingly incomplete. Silicon and the carbon polytypes 
diamond and graphite are particularly relevant because 
long spin relaxation times can be expected in materials 
with inversion symmetry and low atomic number Z. For 
those, the main spin-relaxation mechanism at high tem- 
peratures is the Elliott- YafetfEY) mechanism mediated 
through spin-orbit coupling [J 0] , which scales as Z 2 . 

Silicon, an attractive potential spintronics material due 
to its compatibility with current technologies, has a rel- 
atively large spin orbit-coupling (44 meV) [3J. Never- 
theless^Lepine's electron spin resonance (ESR) measure- 
ments [J, l5( , recently confirmed for low temperatures by 
Appelbaum et al. @, Q, found for its spin relaxation 
time T\ a value of 7 ns at room temperature, which puts 
it well within the usable range. The experimental sit- 
uation is less clear for carbon, whose lower Z-number 
promises longer spin relaxation times. Diamond with a 
spin-orbit coupling of 13 meV Q is especially expected to 
have a long electronic spin relaxation time, which how- 
ever has never been measured. For graphite, the most re- 
cent experimental data from 1961 [9j, [l(| suggest a large 
range for T\ between 1-300 ns. For both diamond and 
graphite, no reliable theoretical predictions have been re- 
ported. Finally, the E~ 2 dependence of the spin-mix am- 
plitude found for III-V semiconductors (lll | is frequently 
assumed to be transferable to other systems such as ele- 
mental semiconductors 12| , whereas proof for that or a 



rigorous calculation of the actual dependence is still lack- 
ing. All this makes truly predictive theoretical methods 
highly desirable to quantify the achievable spin relaxation 
times for delocalized electrons in group-IV materials. 

In this letter, we present such a method, which is 
able to calculate spin relaxation times fully from first- 
principles without adjustable parameters and without re- 
strictions concerning the nature of the band gap. This 
method represents a new reliable and unbiased way to 
calculate spin relaxation times where eigenstates and 
phonon dispersions are calculated self-consistently and 
especially adds to the previously existing body of work 
the capability of calculating parameter-free the effect of 
arbitrary impurities or defects on spin-transport. 

We focus on the EY mechanism, which dominates the 
technologically relevant temperature range above 100 K 
[5|. In the EY mechanism, spin relaxation is induced by 
momentum scattering off impurities or phonons. Within 
the Born approximation, the EY spin relaxation time can 
be related to the momentum relaxation time (which is 
proportional to the carrier mobility) (l3| . The underly- 
ing theory to connect them exists on a phenomenological 
level for III-V semiconductors with direct gap [13. Il5|. 
but not within a first principles framework that includes 
indirect band gap semiconductors 11 1. 



Since a methodology based on density-functional the- 
ory (DFT) to calculate electron mobilities (and thus mo- 
mentum relaxation times) has been recently developed by 
one of us [l6| , what is left to show here is establishing the 
relationship between the spin-flip and momentum scat- 
tering matrix elements. Other than most previous work, 
we do not employ a semiempirical k-p representation of 
the band structure to model the effect of spin-orbit cou- 
pling on the electronic wave functions [13] , but rather use 
the spin-dependent DFT wave functions directly. 

In the presence of spin orbit coupling, Bloch states are 
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given by a mixture of spin- up and spin-down states [l[ , 

*k„ t (r)= [ak»(r)| t> + &kn(r)| I)] e' kr 
*k„4W= K kn (r)| |) - 6*_ k „(r)| f>] e ik ' r (1) 

with lattice momentum k, band index n, and effective 
spins or pseudospins -ft and JL Using Eq. ([I]) for potentials 
that are slowly varying in space on the scale of the unit 
cell m [lH , the relationship between the spin- flip matrix 
elements and the momentum matrix elements becomes 
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with a kj „(G) and & k>n (G) being the Fourier transforms 
of a k , n (r) and 6 k ,n(r) from Eq. ([l}. is a scattering op- 
erator which can refer here to either electron-phonon (i.e. 
lattice) or impurity scattering. Explicit expressions for 
these two types of scattering mechanisms will be given 
below. We assume there is no spin imbalance between 
up and down electrons. The momentum relaxation times 
can be calculated using DFT and density functional per- 
turbation theory (DFPT) [l|| as shown in 0. The 
flkn(G) and b^ n (G) coefficients are obtained from a DFT 
calculation that includes spin-orbit coupling. 

The average spin relaxation time T\ is given by 0, E| 
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where e = e k ,n are the energy bands, Sf the Fermi en- 
ergy, and /(ffk.n) the Fermi function. For phonon scat- 
tering, Ti(k, n) is given by Q 
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with phonon energies w qA and Bose-Einstein occupation 
factors n qA . The temperature dependence of T\ is given 
by the latter two functions, q and A are the phonon 
wavevector and polarization. The spin-flip electron- 
phonon coupling function g is given by [ljj 
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with atom mass M , phonon polarization vectors e qA and 
atom displacements u qA from their equilibrium positions. 
The electron-phonon interaction potential V e - P is given 
by Eq. 19.1 of ref. 

V e _ p (r, p, cr, R + u qA ) = V(r, R + u qA ) 

r V r V(r,R + u qA ) x p- a (6) 
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where the first term is the screened one-electron potential 
V that depends on the electronic position r, equilibrium 
atomic position R and displacement w qA . The second 
term is the spin-orbit potential (Vgo)- p is the electron 
momentum and a are the spin Pauli matrices. Force con- 
stants and phonon fre quen cies are computed from first 
principles using DFPT [H]. 

The inverse spin relaxation time for a density of 
defects or impurities is given by 
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where V is the volume of the supercell. Within the Born 
approximation the spin-flip scattering matrix element is 
given by M„„-(k ft,k' J|) = (nk ff \AV\n'k' J|), where 
the self-consistent scattering potential AV is the differ- 
ence between the potentials of a system with a defect or 
impurity and a reference "unperturbed" system. Elec- 
tronic screening and computational issues due to the use 
of supercells have been addressed as in [lij]. All neces- 
sary DFT and DFPT calculations were performed within 
the local density approximation in Quantum-ESPRESSO 
[20| , using norm-conserving pseudopotentials that explic- 
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itly include spin-orbit coupling [2l|, l22|. We have used 
40x40x40 and 48x48x48 k-point meshes for silicon and 
diamond, respectively and 48x48x16 for graphite, which 
was the maximum possible with our computational setup. 
Convergence was checked by using lower meshes. We find 
that for silicon the convergence error is less than 1 %, for 
diamond 5 %, and for graphite 10 %. 

Since the spin-orbit coupling enters the equations di- 
rectly through the wave functions and not through a 
combination of band gap and spin-orbit splitting param- 
eters, it is a priori unclear to what extent the known 
local density approximation (LDA) band-gap problem is 
a problem for the accuracy of our parameter-free method 
(LDA predicts band gaps significantly smaller than ex- 
periment 23|). Our results show that, for the cases where 
we can compare to experiments, even having a smaller 
LDA band gap leads to good agreement with experimen- 
tal data. By applying pressure we are able to theoret- 
ically modify the direct and indirect band gaps in sil- 
icon. The ab-initio pressure coefficients thus obtained 
(dEg/dP = —1.5 and 0.3 meV/kbar for indirect and di- 
rect band gaps, respectively) are in very good agreement 
with experiments |24| . Artificially changing the band gap 
by varying the lattice strain led to a slow dependence 
(Eg 67 ) of the spin flip amplitude as a function of abso- 
lute gap in silicon. This result gives the opposite trend 
and its dependence is much weaker than the predicted 
scaling of E~ 2 given by perturbation theory applied to 
direct-bandgap semiconductors [Uj, which in the litera- 
ture has been assumed to be correct for Si as e.g. stated 
in Ref. 121. In contrast, we found in the case of diamond 



a dependence of E~ 2 25 . The opposite trend of band gap 
as a function of pressure in the case of diamond is re- 
lated to the absence of d states with the same quantum 
number as the low lying s and p orbitals in carbon [25| . 
In silicon the d states mix with the s and p states thus 
lowering the energy of the conduction band near the X 
point. With increased pressure, these mixed states go 
down in energy and thus give a negative pressure coef- 
ficient. This mixing does not happen in diamond. As 
a result, the pressure coefficient near the bottom of the 
conduction band is positive. 

A calculation of Ti for silicon due to phonon scattering 
using empirical pseudopotential and bond charge models 
has been reported recently by Cheng et al. [12(. Cheng's 
results agree very well with exp eriments. Additionally, 
two very recent papers (26l . |27j used analytical models 
to describe the symmetry of the electron spin-phonon in- 
teractions in silicon in detail. In contrast to that, the 
use of fully first-principles DFT and DFPT here allows 
truly predictive parameter-free calculations and addition- 
ally enables calculation of impurity scattering effects, for 
which no previous work exists. 

Silicon— Our results for the separate contributions to 
the spin relaxation time from acoustic and optical phonon 
scattering are shown in Fig. DJa). Acoustic scatter- 
ing is the most relevant phonon scattering mechanism 
throughout most of the temperature range. In contrast 
to previous suggestions 0, we however find that optical 
scattering starts to become competitive near room tem- 
perature. A T temperature dependence is found for 
phonon scattering. We find that the inclusion of the spin- 
orbit term (Yafet term) in the electron-phonon coupling 
function (second term in Eq. ([6])) is crucial for obtaining 
the right temperature dependence as was first pointed 
out by Yafet Q. Indeed, the cancellation of the zero 
and first order terms between the gradient of the electro- 
static potential and the gradient of the spin-orbit poten- 
tial leads to a higher Ti and much stronger temperature 
dependence. Without the Yafet term (black curve in Fig. 
[TJa)), T\ follows a 7 1 " - 35 dependence which leads T\ to 
be ^3500 times smaller at 80 K and ~30 times smaller at 
room temperature. The calculated temperature depen- 
dence of T~ 3 - 5 lies between the high temperature limit 
of T~ 2 5 and low temperature limit of T -5 calculated by 
Yafet for acoustic phonons. Also in Fig. HJa) we show 
the results for spin relaxation times for impurity scat- 
tering as a function of temperature for different carrier 
concentrations. Increasing the carrier concentration by 
adding more donors in the examined range makes the im- 
purity scattering mechanism to become relevant at tem- 
peratures below 150 K, lowering the total spin relaxation 
time as expected. 

We have used Matthiessen's rule to add the lattice and 
impurity contributions to the total spin relaxation time. 
In Fig. [Ub) we compare our results with the ESR exper- 
iments of Lepine and Appclbaum's recent measurements. 
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FIG. 1. (Color online) (a) Phonon contribution to the elec- 
tronic spin relaxation time Ti in Si for all (red squares), 
acoustic (blue diamonds) and optical (green squares) phonons 
with Yafet term, and all phonons without Yafet term (black 
circles). Dashed curves are contributions to Ti from P - 
impurity scattering at carrier concentrations of 7.4 x 10 14 (+), 
7.8 xlO 15 (x) and 8xl0 16 cm -3 (*). (b) Theoretical Ti (lines) 
for Si compared with experimental values [fj for carrier con- 
centrations of 7.4 x 10 14 (circles), 7.8 x 10 15 (squares), and 
8 x 10 16 cm -3 (diamonds) and with [?| (triangles). 



We find very good agreement in the region above 150 K, 
where the EY mechanism is considered to be dominant. 
At room temperature the calculated relaxation time of 
T\ = 4.3 ns is well in the usable range (l2j . 

To check if our method would also produce useful re- 
sults in the degenerate-doping regime, we also calculated 
the room-temperature spin relaxation time for a donor 
concentration of 1.8 x 10 19 cm~ 3 , for which Dash et al. 
recently suggested a lower bound for the electronic spin 
lifetime of initially 140 ps [H], then 285 ps 0, nearly 
two orders of magnitude smaller than Lepine's results in 
the low-doped regime. This indicates the relevance of 
impurity scattering in lowering the spin lifetime in the 
highly-doped regime, which is important to keep deple- 
tion zones reasonably small. We find a value of 180 ps, in 
the range of the experimental values. From this compar- 
ison, it seems that the changes in band structure (impu- 
rity bands) and wave functions from degenerate doping 
[30j play a lesser role for the value of the spin relaxation 
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FIG. 2. (Color online) Spin relaxation due to phonon scatter- 
ing in graphite. Red squares (black circles) were calculated 
with (without) inclusion of the Yafet term. 



time, since they are not included in our theory. We also 
find a weak dependence of the spin-mix amplitude as a 
function of band gap. 

Graphite— The metallic character of graphite requires 
integration with special attention to the Fermi line that 
goe s from the K to the H point of the Brillouin zone 
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and inclusion of a minimum of two bands into the 
calculations. We obtain a spin relaxation time of ~5.8 
ns at 300 K (Fig. [2]). This result is at the lower end 
of the available experimental data 0, [Io| which report 
values between 1 and 300 ns. We find T\ to be strongly 
anisotropic. By just considering scattering along the c 
axis, we get a much lower T\ of 0.1 ns. We obtain a strong 
temperature dependence of T -4 - 5 , close to the theoretical 
low-temperature limit Q , as expected from the large in- 
plane Debye temperature of 2500 K. Without the Yafet 
term, T\ would be ~110 times smaller at room tempera- 
ture and follow a linear temperature dependence. 

Diamond— Figure [3] shows the calculated T\ for dia- 
mond considering only phonon scattering. As expected, 
diamond has a considerably larger T\ than silicon with a 
value of 180 ns at room temperature. However, it is sig- 
nificantly smaller than the 10 2 — 10 4 s speculated recently 
[32j using a formula [33| that in the case of silicon gives 
T\ = 1 s at room temperature, eight orders of magnitude 
larger than experiment, as had already been cautioned in 
Ref. 33( . This formula is based on the Waller theory of 



relaxation by modulation of dipole-dipole interactions by 
lattice vibrations 3^ . This mechanism is very weak and 
typically results in spin relaxation times many orders of 
magnitude larger than measured experimentally [33j . On 
the experimental side, data are available for the dephas- 
ing times for the well-studied nitrogen- vacancy center in 
diamond (us to ms) [H, [35| , but an experimental study 
of the spin relaxation of conduction electrons is still lack- 
ing. In order to check if the magnitude of our results is 
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FIG. 3. (Color online) Calculated phonon contribution to 
the electronic spin relaxation time in diamond from all 
(red squares), acoustic (blue diamonds) and optical (green 
squares) phonons. The black circles show spin lifetimes cal- 
culated without the Yafet term. 



reasonable, we have performed mobility calculations [36[ 
of conduction electrons at room temperature obtaining a 
value of 130 cm 2 /Vs, which is within the range of recent 
experimental Hall data H3 (100-660 cm 2 /Vs). Since T\ 
for Si agrees well with experimental data and the electron 
mobility for Si has been accurately predicted previously 
fl(| . the link between the matrix elements that deliver 
the mobility and the spin- relaxation time (Eq. @) has 
been validated. Based on this chain of benchmarks, we 
arc confident that the predictions for T\ in diamond and 
graphite are sensible, and will also be reliable for other 
materials where spin relaxation is dominated by the same 
processes. As temperature decreases, the T~ 5 temper- 
ature dependence in diamond (expected as in graphite 
due to the high Debye temperature of 2200 K) leads to 
a much larger T± than in silicon (^110 times larger at 
150 K). Accordingly, scattering from acoustic phonons 
is dominant throughout the whole temperature range up 
to room temperature. When we do not include the Yafet 
term in the calculation, T\ becomes ~2300 times smaller 
at 300 K and with almost no temperature dependence 
(~ illustrating again the relevance of including 

this term in the electron-phonon coupling function. Since 
the calculated band gap dependence of the spin-mix am- 
plitude is Eg 2 - 25 and the LDA band gap is 20 % smaller 
than the experimental gap, we estimate the true T± to 
vary within a factor of 1-1.6 from the reported ab-initio 
values. 

In conclusion, a new parameter-free first-principles 
method to obtain spin relaxation times for phonon and 
impurity scattering has been presented, which is gener- 
ally applicable to arbitrary systems. We have bench- 
marked this method for the Elliott- Yafet dominated tem- 
perature regime in silicon with very good agreement with 
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experiment for phonon and impurity scattering and a 
limit for T\ of 4.3 ns at room-temperature. For silicon, 
we find a weak dependence of the spin-mix amplitude as 
a function of band gap. We predict a stronger T~ 5 tem- 
perature dependence for phonon scattering in diamond 
with T± = 180 ns at room temperature. This value is 
significantly smaller than the 10 2 — 10 4 s recently esti- 
mated [H] using a formula [HI that finds equally for Si 
a value that is eight orders of magnitude higher than ex- 
periment. Although no experimental results exist for Ti, 
electron mobility measurements support our results. For 
graphite we find a T -4 5 temperature dependence with a 
T\ of 5.8 ns, at the lower end of the experimental range. 
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